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Abstract 



We study a multi-species exclusion process with inhomogeneous hopping rates. This model is equivalent 
to a Markov chain on the symmetric group that corresponds to a random walk in the affine braid arrangement 
[11] . We find a matrix product representation for the stationary state of this model. We also show that it 
is equivalent to a graphical construction proposed by Ayyer and Linusson 4 , which generalizes Ferrari and 
Martin's construction [8]. 

1 Introduction 

The coupling of randomness with algebraic or arithmetic structures can lead to beautiful combinatorial results. 
The study of random partitions of integers and its extension to three-dimensional plane partitions [T7l HJ [15] 
provides a prominent example. Recently, Lam studied Markov chains that can be represented geometrically as 
random walks on a regular tessellation of a vector space on which an affine Weyl group acts [TT] . When this Weyl 
group is a symmetric group, the stationary distribution of the chain displays remarkable combinatorial properties 
that were further explored by Lam and Williams, leading to various conjectures |12j . The Markov chain studied 
in |12j is equivalent to a multi-species exclusion process with inhomogeneous transition rates as shown by Ayyer 
and Linusson in a very recent publication [5] . This suggests that the powerful techniques that were developed in 
non-equilibrium statistical mechanics to analyze the asymmetric exclusion process [14, 6, 5 should be relevant 
to this more mathematical problem. Indeed, Ayyer and Linusson conjectured that the stationary state can 
be obtained by generalizing the elegant graphical algorithm, invented by Ferrari and Martin [5] to solve the 
homogeneous iV-species totally asymmetric simple exclusion process (N- TASEP); in [4], they generalized Ferrari 
and Martin's algorithm and explored the consequences for Lam and Williams' conjectures. 

In the present work, we solve the TV- TASEP with inhomogeneous transition rates by using a generalized 
matrix product Ansatz. More precisely, we show that a suitable deformation of the homogeneous algebra studied 
in [2] allows us to calculate the stationary state of the inhomogeneous N- TASEP. We also show that our 
matrix product solution is equivalent to the algorithm conjectured in [4]. 

The outline of this work is as follows. In section [2] the model is defined. In section [3] we give a matrix 
product solution to the stationary state. In section |4j we explain the graphical construction of the stationary 
state conjectured by Ayyer and Linusson, and we show that this construction is equivalent to the matrix product 
solution. We give concluding remarks in section [5] 

2 Definition of the model 

We consider an L-site periodic chain in which each site takes a non-negative integer value 1,2,..., or N + 1 
(this is called the N-species problem). Each pair of nearest neighbor sites exchanges their values according to 
the following continuous time stochastic dynamics: 



a b — > b a 




(a < b) 
(otherwise) 



(2.1) 
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The special case x a = 1 for all a corresponds to the homogeneous A-TASEP [2 |3 [H [7] . When the transition 
rates a; a 's are not equal to each other, we say that this dynamics is inhomogeneous (although one could consider 
even more complicated or general transition rules). When L = N + 1, and there is exactly one particle of each 
species in the system, the process is equivalent to Lam and Williams' process on the symmetric group Sl [H]- 
The inhomogeneous A-TASEP is governed by the master equation 

||P)=MW|P) (2.2) 

for the vector \P) = ^P(ji • ■ • • ■ ■ Jl), where P(ji ■ ■ ■ jt) is the probability of finding the system in 

a configuration ji ■■■jj J . The generator matrix (Markov matrix) M^ N ^ is the summation of local operators 

( A^Loc ) that acts on the spaces corresponding to ith and (i + l)st sites of the chain: 

L N+1 , 

MW=^(MW) , <J= e(a-b)(\ba){ab\-\ab}{ab\) with 6 (a - b) = (2.3) 

i=l 1 a,b=l ' L V - 

We write m = (mi, . . . , mjv+i) (with m, G Z>o) for the sector that contains particles of type k. Because 
of the conservation of the number of particles of each type, we have the decomposition = M m . In 

particular, we consider basic sectors, i.e. m, > 0. Since the identification of local states N + 1 — > N maps the 
A-species dynamics to the (N — l)-species dynamics, we have the spectral inclusion: the spectrum of the sector 
m = (m 1; . . . , 771/v, TOjv+i) contains that of the sector m! — (mi, . . . , m^v + m w+i)[^] This inclusion relation 
indicates that there exists a "conjugation matrix" ip m such that 

M m ^ m = ^ m M m ,. (2.4) 

This matrix allows one to lift states that belong to the (A — l)-species sector m! and to construct eigenstates 
for the A-species sector m. In particular, the stationary state of the A-TASEP (the kernel of AfW) can be 
constructed recursively if one knows the sequence of the conjugation matrices. 



3 Generalized matrix product Ansatz 

The idea of the matrix product Ansatz is to express the probability of finding each configuration ji ■ ■ ■ jh as a 
trace over a suitable algebra 

^••■i^^K'-ljf) (3-D 

with a normalization constant Z. The operators Xj's must satisfy suitable algebraic relations, which are 

usually infinite dimensional [51 [S] . For the homogeneous A-TASEP, the operators X { ^ ] were constructed recur- 
sively [7J [THJ as tensor products of four fundamental operators 5, e, A and 11 that act on an infinite dimensional 
space A = ® At >o C|yu)) as 

= { t>t 4m>>Hm + i»,4m>> = {P £ = $ %» = lM»- (3.2) 

One can easily verify that the following quadratic relations are satisfied: 

Se =l, 5 A = 0, Ae = 0. (3.3) 

In [16], the stationary state for the homogeneous A-TASEP has been expressed as a matrix product form, using 
an algebraic interpretation of Ferrari and Martin's algorithm [8]. More recently the matrix product Ansatz 

1 More general inclusion relations are satisfied for the homogeneous case, see [3]. 



2 



technique was generalized to obtain a conjugation matrix that satisfies "conjugation relation" (2.4) between 
systems having different numbers of species [2]. We now explain how the same ideas can be adapted to the 
inhomogeneous case. 

The stationary state for N = 1 case is trivial, i.e. all the possible states in each sector are realized with a 
same probability. This can be regarded as one dimensional representation of the matrices — = 1. 

3.1 The 2 species case 

For N — 2, a matrix product stationary representation has been known, even in the inhomogeneous case [S]. 
One possible choice for the algebra is 



Y {2) y (2) _ 1 y(2) 



X 



(2) 



3 > 



Y (2) Y (2) _ 



X2 



2 ) 



y(2) y(2) _ y(2) 



(3.4) 



Defining an operator-valued vector X^ 2 ) = 



2) 



1 





(2) 
2 

(2) 



we observe that the following decomposition exists 



5 
A 



y x e y 1 B + y 2 A 




(3.5) 



where we have set y\ — l/x\,y2 — 1/2:2 and the matrices <5, e, A satisfy (3.3 ). This decomposition can be written 
more formally as the product of two rectangular operator-valued matrices and a^ 2 ) of sizes 2x1 and 3x2 
respectively 



X (2) 




(3.6) 



where the symbol * means that we perform tensor products amongst the elements of the matrices a^ 2 ) and 
sS x \ More generally, for matrix-valued matrices Y = (Y uv ) uv , Z = (Z uv ) uv , it represents the product Y * Z = 
(%2 W Y UW ®Z WV ) . Since and a^i are scalars here, the tensor product reduces to the ordinary product. 



where 



we observe that the matrices a^ 


and 


a^ 2 ) satisfies 








a«< 


= 


a«$ 


Da« - 


M^W 2 ) - 


a< 2 )( 




8< a >€ 


)a( 2 ) - 




f 


), s (2) = 


(• 


i) 



i«®a« 

,( 2 W 2 ) 



(3.7) 
(3.8) 



(3.9) 



Since = 0, equation (3.7) is the relation in the usual matrix product Ansatz [5]. The relation (3.8 1 implies 



that the matrix i() m whose elements are given as 

O'l ■ • -jL\lpm\kl --'kt) 



Tr 



( (2) 



(2) 



(3.10) 



intertwines the dynamics of the sectors m = (mi, W2, TO3) and m! = (mi, TO2+TO3), i.e. it satisfies equation (2.4) 
Formally we write the stationary state for the sector m! as ipm' with (ji ■ ■ ■ :/l|Vw|1 •■■!)= Tr^a^ ■ ■ • aj* iT = 
1. Then the stationary state for the sector m can be rewritten as "0?nVw- 



3 



For a simple nontrivial example, the stationary state of the sector (1, 1, 2) is given as 



^(1,1,2)^(1,3) 





1222 


2122 


2212 


2221 




1233 


/ vl 











\ 


1323 


vl 


2/12/2 










1332 


vl 


2/12/2 


vl 







2133 





vl 


2/12/2 


vl 




2313 








vl 


2/12/2 




2331 











vl 




3123 





vl 










3132 





III 


2/12/2 







3213 


vl 





vl 


2/12/2 




3231 


2/12/2 








vl 




3312 








vl 







3321 


\ 2/12/2 


vl 





2/1 


/ 



/ vl \ 

2/2 (2/1 +2/2) 



/ 1 \ 
1 

1 

V 1 / 



vl + 2/12/2 + 2/1 



2/12/2 
2/2 (2/1 + 2/2) 
vl 
vl 



vl 



vl 



(3.11) 



2/2 (2/1 + 2/2 

2/12/2 + 
2/2 (2/1 + 2/2) 
2/1 

V 2/i + 2/12/2 + vi } 
3.2 The iV-species case 

We now explain the generalized matrix product Ansatz for general values of N (see [5] for more details). 
Reversing the construction in the last subsection, we start with the following relation that we call "hat relation" : 

M^a^WW - a^W^M^ = S (JV) ®aW - a^Sa^ , (3.12) 
where a^ and a^ are operator-valued matrices of size (N + 1) x TV. We write their elements as (j\a^ N '\k) = 
CLjk \ (j\sS N ^\k) — otj^K We know that, if we can construct a couple {aW,a^} (for the general integer of N) 
that satisfies (3.12), the matrix tp m defined as 

Vl ■ ■ ■ jL^mlki ■■■k L )= Tr( a W • • • afX) (3-13) 

satisfies the conjugation relation (2.4). Here the configurations j\---jh and k\---kh belong to the sectors 
m = (mi, . . . ,rnjv+i) and m! = (mi, . . . ,rriN + mjy+i), respectively. Furthermore, the stationary state of the 
sector m can be written by the product of conjugation matrices 

\P)m = VVnVw ' ' ' 1p(mi,L- mi ) (3-14) 

if all of them are nonzero (we note that the conjugation matrix ip m lifts up other eigenstates from lower sectors 
as well). 

We set iji — 1/xi and define the operator B = 1 — A. An explicit solution to the hat relation (3.12 1 is given 

by 



AN) _ 

l jk 



= < 



A®(j-V®5®I®( k - j - 1) ®e®l® (N - k - 1 '> 

(l£=i V i A 9 (*- 1 '>®B®l 9 <- k - i - 1) + y k A® {k -^ 
Eti" 1 y l A^-V®B®l® ( - N - l ~V + Vn a®( n -V 


l®(>°-i)® e (g l l®(N-k-i) (fc < j — l — TV) 

l^- 1 ) (k = j — 1 = N) 

(otherwise). 

(N) 



(JV-fc-l) 



(j<k< TV), 
(j<k = N), 

U = k), 

(k<j-l = N), 

(j-l = k = TV), 
(otherwise) 



(3.15) 



(3.16) 



The difference from the homogeneous case appears in a) fc for j = N + 1. Indeed we retrieve the solution to the 
homogeneous case [HE] by setting yj = 1. We change the definitions of a^ and a^ as a^ = 



1 

Vi 



a« = 



for compatibility with the general forms (3.15), (3.16) 
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By a direct calculation, one can show that equations (3.15), (3.16) give a representation to the algebra 



defined by the hat relation (3.12) i.e. 



XjQjkQ/ji k' Xki^Cljk'Clj'k 


— ajkClj'k') 


— CljkClj 


k' 


Clj kClj 


k' 


(j<j'Ak< k'), 


—Xk(a,jk>a,jik 


— ajkdj'k') 


— CljkClj 


k' 


CljkClj 


k' 


(j = j' Ak<k'), 


^Cj'Clj'kCljk 1 -^k{cijk'Clj'k 


— ajkdj'k' ) 


— CljkClj 


k' 


— cijkcij 


k' 


(j>j'Ak< k'), 




~% jCtjk&j' k' 


= ajkcij 


k' 


— cijkaj 


k' 


(j<j'Ak>k'), 







= ajkCbj 


k' 


— cijkaj 


k' 


U=j' Ak>k'), 






CljkClj 


k' 


Cijkaj 


k' 


(j>j'Ak> k'). 



(3.17) 
(3.18) 
(3.19) 
(3.20) 
(3.21) 
(3.22) 



It is straightforward to prove that the above relations are satisfied simply by substituting equations (3.15) and 
(3.16). In appendix, we prove the first identity (3.17) as an example. For N — 3 and 4, the result can be written 



explicitly as 



,(3) _ 



/ 1®1 



2/1 e 



1 



<5®e 




V 



5® 1 

A® 5 

A®A 
y\B®t 
+y 2 A®B 
+y 3 A®A J 



( 



a< 3 > = 









e(g)l l®e 



\ 



/ i®i<g>i 
o 




yie<g>l<gl 








e<gl<g)l 



<5<g>e(g)l 
A®\®\ 





y 1 B<Z>e<S)l 
+y 2 A®e®l 



5®l®e 
A<g>5®e 
A<S>A<E)t 


y 1 B®l®e 
+y 2 A®B® 
+y 3 A<E>A<8) 








6® 1(g)! 

A®5®1 

A®A®8 

A®A®A 
yxB®!®! 
+y 2 A®B®t 
+y 3 A®A®B 
+y 4 A(E)A<E)A J 







1(E) 1 J 
\ 



(3.23) 



(3.24) 



1; 



(3.25) 



)e®l 1® l®e l(g)l(g)l J 



We now generalize the form (3.5) or (3.6). The form (3.14) can be written as the matrix product form (3.1) 
with the matrices 



(3.26) 



thanks to the "sector specificity"^] 



Let us consider the element (3.13) for configurations i\---jh an d ki ■ ■ ■ k^ of the sectors (mj, . . . , mjy+i) 
and (mi, . . . , mjv + tojv+i) with ji < ji+i- Since cijkaj'j' = for j < j' < N and k ^ k', we need to set ji = ki 
for i < L — tojv+i so that (3.13) is nonzero, and we have 

Tr (an 



■a Y \ci 22 ■ 



a 22 ■ ■ ■ a NN ■ ■ ■ a NN a( N+1 } N ■ ■ ■ a^+i)^) 

N ^nire+lH hnijv+l 



m N+1 
UN 



(3.27) 



Therefore the stationary weight of the configuration ji ■ ■ ■ jl is ]~Iri=i V" h '" JV+1 . As far as we treat basic 

sectors, the matrix product ■ ■ • a j^l L contains = A®( N ~ 1 ' > or 0. This implies that the matrix product 
is transforme d in to a tensor product of the form e c A5 c = \c))((c\ multiplied by a monomial of y n 's (by applying 
(O) and A 2 = A), otherwise it is 0. Since Tre c AS C 



the relation 
of y n 's or 0. 



1, the trace of a 



(AO 



AN) 



jiki 

Note that the trace is not always finite if we consider a non basic sector. 



is a monomial 



2 When fci • •• hi, does not belong to the sector m' , the trace j3 ■ 1 3[ ) is always 0. This specificity is because the numbers of 8't 
and e's are different in the matrix product, see [2] for details. 
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4 Graphical construction of the stationary state 



We have found a solution to the generalized matrix product Ansatz that constructs conjugation matrices of 
the inhomogeneous 7V-TASEP. In this section, we first review the algorithm that A. Ayyer and S. Linusson 
devised [3] to calculate the stationary weights by defining an inhomogeneous extension of the seminal algorithm 
by Ferrari and Martin^. Then we show that the solution to the matrix product Ansatz of the previous section 
is equivalent to Ayyer and Linusson's construction. 



4.1 Ayyer and Linusson's algorithm 

The algorithm of Ayyer and Linusson [4] constructs the stationary state of the A^-species sector (mi, . . . , mjv+i) 
from the (N — l)-species sector (mi, . . . , m N + m N+ i). It is provided by two maps F, W from an (N— l)-species 
configuration and a configuration consisting of black and white boxes, to an ./V-species configuration and a 



polynomial in j/j's. Figure 4.1 is helpful to understand the algorithm. 



(i) Let us set two lines. On the upper line, there are mi + ■ • • + m^r black boxes I and mj\r+i white boxes □ 
as Ci • • • cl (ci = M, □). On the lower line, we give a configuration fci • • • Ul of the (N — l)-species sector 
(mi, . . . , tojv + mjv+i). Thus, on the lower line, there are m„ i/'s (1 < v < N — 1) and (m^v + mjv+i) AT's. 

(ii-1) Let {ii, ■ ■ ■ , £J} be the positions of the mi l's on the lower line. For the first 1 located at find 
the nearest black box Cj< = ■ with i' < and put 1 on it. If there is no such black box, put 1 on 
the rightmost black box. For the second 1 located at i 2 , find the nearest unoccupied black box c,< = ■ 
with i' < i 2 and put 1 on it. If there is no such black box, put 1 on the rightmost unoccupied black 
box. We draw an arrow from ij 1 ' to the targeted black box. (In the following procedure till (h-^), we 



always draw an arrow in the same way, see figure (4.1 1 ) We iterate this procedure mi times, i.e. find the 



nearest unoccupied black box c.j< = I with i' < i^ 1 ' for the £th 1 located at iy\ and put 1 on it or on the 
rightmost unoccupied black box if i' does not exist. 

(2j) (2} 

(ii-2) Let {i\ , . . . , i m2 } be the positions of the m2 2's on the lower line. There are {m-i + ■ ■ ■ + mjv) unoccupied 
black boxes on the upper line. We iterate the following procedure mi times: find the nearest unoccupied 

(2) 

black box c,-/ = ■ with i <i\ for the £th 2 (1 < £ < 7712), and put 2 on it or on the rightmost unoccupied 
black box if i' does not exist. 

(ii-i/) In the same way, we go on for v = 3, 4, . . . , (N — 1). Let (4 , ■ ■ ■ , it!} be the positions of the m„ ^'s 
on the lower line. There are (m v + • • • + mjv) unoccupied black boxes remaining on the upper line. We 
iterate the following procedure m v times: find the nearest unoccupied black box Ci' — ■ with i' < for 
the iih v (1 < £ < m„), and put v on it or on the rightmost unoccupied black box if i' does not exist. 

(iii) There are mjv unoccupied black boxes remaining. Put A^s on them. 

(iv) Put (N + l)s on the m^r+i white boxes. We have thus constructed a configuration F(ci ■ ■ ■ cl, fei • • • /cl) 
of the A^-TASEP on the upper line, belonging to the sector m. 

(v) Define a vector from the stationary state \P m ') of the sector m! as follows: 

\P m ) = W (3i •••3L,ki-~ k L )\ji ■ ■ -j L )(ki ■ ■ ■ k L \P m ,). (4.1) 

The summation runs over ji ■ ■ ■ jl and hi" -k^ belonging to the sectors m and m'. If there exits a 
configuration ci • • • cl such that 

F(c 1 -"C L ,k 1 -"k L )=j 1 -"j L , (4.2) 
the coefficient W(ji ■ ■ -jl, hi ■ ■ ■ k^) is defined as the product of the following weights {wi, • • • , rat}. Draw 



a vertical line between each bond between site i — 1 and i, see figure 4.1 When an arrow connecting two 



G 



(ii-1) 



443123444 
(ii-3) 



t t 






443: 




1 

I 3 



L 



(ii-2) 



■ □■□□■□I 



44312344 4 



(iii)(i 



3 2 



1 2 



5 4 



weights 1 1 y 2 1 y 2 y 3 1 y 4 1 



5 3 



L 



Figure 1: An example of Ayyer and Linusson's algorithm. 



v's on the upper and lower lines in our figures, we say "the value of the arrow is u" : t— | . Each weight Wi 
is defined as 

1 (ji < N), 

ye (ji = N + 1, and I is the minimal value of arrows that cross the ith vertical line), (4.3) 
UN {ji = N + 1, and no arrow crosses the ith vertical line). 

Then we have W(ji ■ ■ • jh,k\- • ■ k^) = rii<i<L w i- If there is no configuration c\ ■ ■ ■ cl such that equation 
(4.2 ) is satisfied, we define W (ji ■ ■ ■ Jl, k\ ■ ■ ■ k£) = 0. 



Ayyer and Linusson conjectured that the form (4.1) gives the stationary state \P m ) of the sector m [4]. 
There the definition of the weight looks different from our W, but is equivalent by multiplying it by a constant. 

Figure 4.1 gives an example of the algorithm for the sectors m = (1, 1, 2, 1, 4), to' = (1,1,2,5), where the 
upper and lower lines are 

ci • • ■ c L = k!---k L = 443123444. (4.4) 

According to the algorithm, the configuration of the sector to and the weight are obtained as 

^(■■□■□□■□■, 443123444) = 325155453, 1^(325155453,443123444) = ylym- (4-5) 



4.2 Equivalence of the matrix representation and Ayyer and Linusson's algorithm 

We explain the relation between the matrix product representation and Ayyer and Linusson's algorithm (follow- 
ing [5] for the homogeneous case). The elements a^' act on a basis vector \fix, . . . , yN-i)) = • • • ®\yN-i)) € 
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A®( N ~V as 



a ) k m 



,Hn-i)) 

Imi> ■ ■ ■ ,H - i 
|Mi> •••iMj - 1 
. . .,Hn-i)) 

2/min{£,/c}lMl) • • • j Mfc 
^|/X1,...,^JV-1» 





, /i fe + 1, . . . ,hn-i)) (j < k < N, hi = ■ ■ ■ = Hj-i = 0,/ij > 0), 
,Mjv-i)) (j < k = N, fii = ■ ■ ■ = fij-i = 0,/ij > 0), 

(j = fc,/xi = • • • = ;U.,_i = 0), 
1, . . . ,/xjv-i)) (& < j - 1 = N,(ii = ••• = M-i =Q,He> 0), 
(j - 1 = k = N, hi = • • • = Ht-i =0,Ht > 0), 
(otherwise) . 



(4.6) 



Let ji • • ■ jh and ki • 
matrix product a,j 1 k 1 



■ k^ belong to basic sectors m and ml '. If the vector \fii, - ■ ■ ,Hn-i)) is not killed by a 
• ■CLj L k L > the series of matrices {aj 1 k 1 , ■ ■ ■ , a j L k L } §i ve a series ("trajectory") as 



|Mi> • 



(JV) 

1 — > v l \hi,. 



■ • ) MiV-l) 



(JV) 

2 ii_i*L_i . £-1 • 

I > V L - 1 V L \H l ,...,fi N -l. 



■ > MjV-l) 



31 fe ,l 



ui • • ■ • ^Miv-i))) 



(4.7) 



where we have set 



(JV) 



fcJMi 



,/$"-!» =«il/4,---,/4-i)) (m£ +1 = A*v). 



(4.8) 



Since a^ fc ' increases H> ^tfe 

= ^} for v < N - 1, we hnd 



1 and decreases //j i-> /Zj — 1 (see the action (4.6)), and we have = v} = 
\hi,---,HN-i))- In ot 



' MiV-l) 



rer words, . . . , hn-i}} is an 



eigenvector of the matrix product with a nonzero eigenvalue. If a trajectory (4.7) is given, one notices that it is 

= for (Mi,-.-,Mat-i) ^ (Ml.' 



j (JV) 

unique and ctj lk \ 



AN) 



I Ml) 



) Miv-i/ 



. , Hn-i)' We regard in the vector 



. . . , fiN—i)) as the number of arrows with value The trajectory gives one graph of arrows which is the 
same as obtained by Ayyer and Linusson's algorithm since the local action (4.6) is compatible with it. (For 
example, when j < k < N , the action of decreases # of arrows with value j and increases # of arrows 
with value k. This is allowed only if there is no arrow with value v < j, otherwise it kills the vector.) We can 
also see a compatibility of the local action for the coefficient Vi = Wi as well as arrows, and thus the nonzero 
eigenvalue is identical to W(ji ■ ■ ■ jl, fci • • • 



Tr (a 



AN) 



AN) 



= W(ji---j L ,ki---k L )- 



(4.9) 



When there is no nonzero eigenvalue, this equation is also true since W(j± ■ ■ • jt, fci •• ■ kij) = 0. 

For example, for configurations ji ■ ■ ■ jt — 325155453, ki ■ ■ ■ k^ = 443123444 as in figure |4.1| the matrix 
product aj 1 k 1 ■ ■ ■ a>j L k L gives a trajectory 



2 (4) 
'31. 



|0,0,1)) ^> |0,0,. 



y 4 |0,0,( 



2/22/32/4|0, 1, 1)> 



|0,0,( 
y 2 2 y 3 y 4 |0,l,2)> 



i (4) 



2/ 3 y 4 |0,0,i)) 



y 2 2 y 3 y4|o,o,2)) 



2/2^32/410,1,1) 

2/ 2 2 2/32/ 4 |0,0,l)) 



(4.10) 



and we have a jlkl ■ ■ ■ a jL k L |0, 0, 1)) = y^Vsy^O, 0, 1)), and we have Tr(a jlkl • • • a jhkL ) = yly-$y\- 



5 Concluding remarks 

In this work, we applyed the generalized matrix product Ansatz to represent the stationary weights of the 
inhomogeneous A^-species TASEP. We also explained that our solution to the Ansatz is equivalent to Ayyer and 
Linusson's combinatorial algorithm. Our analysis was motivated by some conjectures proposed by Lam and 
Williams 12 . For example, for L = N + 1, they claim that the stationary probability is a polynomial with 
respect to the hoping rates with non-negative integer coefficients and is a non-negative integral sum of Schubert 
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polynomials. The first observation is an obvious consequence of Kirchhoff 's matrix tree formula |18j . (Note that 
this is also an outcome of our matrix product solution because all operators have positive entries.) However, the 
relation with Schubert polynomials is still unclear. Another interesting problem would be to extend the present 
study to the partially asymmetric case and to classify the multi-species systems with arbitrary inhomogeneous 
hopping rates that can be solved by the matrix product Ansatz. Finally, we note here that the inhomogeneities 
are linked to the particles rather than to the underlying lattice. Extending our approach to models with lattice 
defects (such as the Janowsky and Lebowitz model in which the insertion of a slow bond can generate a shock 
[9]) remains a very challenging open question. 
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A Proof of equation (3.17) 



Here we show the first case ( |3.17 ) of the algebra. 
The case when j ^ k, f < N . We have 



dj K = or jth component of a,j K is S (for k = k,k'), 
cij'K — or jth component of aj> K is A (for k = k, k'). 



(A.l) 
(A.2) 



In any cases, we find ajkdj'k' = aj'kO>jk' — 0, and thus the left-hand side is 0. The right-hand side is also 
thanks to ajk = cij>k' = 0. 



The case when j = k, j' < N . The left hand side is —Xjdjk'a-j'j — thanks to ayj = 0. The right-hand 



side is again thanks to ajk — aj'k' = 0. 
The case when j < k, j' = N + 1 . Since 



a jk a fk , = a jk ,a fk = y j A®( j -V®6m® ik - j - 1) ®em® ik '- k - 1) ®e®l® {N - k '- 1) , 
the left hand side is calculated as 

- A^ j - 1 ^8^ (k - j - 1) ®em 0(k '- k - 1) ®e^ N - k '- 1 \ 
which agrees with the right-hand side. (We read • • • l^-^igie®!^" 1 ) = • • • l^-k- 1 ) for k' = N.) 

. The left-hand side is calculated as 



(A.3) 



(A.4) 



The case when j = k, j' = N + 1 

- rj,, , = _A®y-l)gl®(*'-i) e ®l®( JV -*'- 1 ) 
which agrees with the right-hand side. 
The case when 



(A.5) 



j > k,j' = N + 1 . We have ajkdj'k' = thanks to ajk — 0. Since the fcth component of 
each term of aj'k is e, and ajk>=0 or the kth component a^y is A, we also have aj k iaj' k = 0. Thus the 
left hand side is 0. The right hand side is also thanks to ajk = fljfe = 0. 
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